Double migration of the endangered Tricyrtis formosana (Liliaceae) in Japan.

The Ryukyu Islands of Japan are a biodiversity hotspot due to geographical and historical factors. Tricyrtis formosana is a perennial herbaceous plant that commonly found in Taiwan. But only a few populations have been identified in a limited habitat on Iriomote Island, while populations of unknown origin occur near human settlements in an area on the main island of Okinawa. To better understand these populations of the phylogenetic uniqueness and intrinsic vulnerability, we conducted comparative analyses including (1) phylogeny and population structure with MIG-seq data, (2) photosynthesis-related traits of plants grown under common conditions and (3) transcriptome analysis to detect deleterious variations. Results revealed that T. formosana was split into two clades by the congeners and that Iriomote and Okinawa populations independently derived from ancestral Taiwanese populations in each clade. Photosynthetic efficiency was lowest in the Iriomote population, followed by Okinawa and Taiwan. Transcriptome analysis showed that the Iriomote population accumulated more deleterious variations, suggesting intrinsic vulnerability. These results indicate that each T. formosana population in Japan is phylogenetically unique and has been independently dispersed from Taiwan, and that the Iriomote population presents a high conservation difficulty with a unique photosynthesis-related characteristic and a larger amount of deleterious variations.

Phylogeographical evidence for historical long-distance dispersal in the flightless stick insect Ramulus mikado.

Exploring how organisms overcome geographical barriers to dispersal is a fundamental question in biology. Passive long-distance dispersal events, although infrequent and unpredictable, have a considerable impact on species range expansions. Despite limited active dispersal capabilities, many stick insect species have vast geographical ranges, indicating that passive long-distance dispersal is vital for their distribution. A potential mode of passive dispersal in stick insects is via the egg stage within avian digestive tracts, as suggested by experimental evidence. However, detecting such events under natural conditions is challenging due to their rarity. Therefore, to indirectly assess the potential of historical avian-mediated dispersal, we examined the population genetic structure of the flightless stick insect Ramulus mikado across Japan, based on a multifaceted molecular approach [cytochrome oxidase subunit I (COI) haplotypes, nuclear simple sequence repeat markers and genome-wide single nucleotide polymorphisms]. Subsequently, we identified unique phylogeographic patterns, including the discovery of identical COI genotypes spanning considerable distances, which substantiates the notion of passive long-distance genotypic dispersal. Overall, all the molecular data revealed the low and mostly non-significant genetic differentiation among populations, with identical or very similar genotypes across distant populations. We propose that long-distance dispersal facilitated by birds is the plausible explanation for the unique phylogeographic pattern observed in this flightless stick insect.

Temporal and interspecific dietary variation in wintering ducks in agricultural landscapes.

Farmlands are becoming more important as waterfowl foraging habitats, while natural wetlands are being lost globally. However, it is unclear how waterfowl coexist in agricultural landscapes by resource partitioning. We evaluated the diets of seven sympatric dabbling ducks foraging in rice paddy and lotus fields around Lake Kasumigaura, the second largest lake in Japan, during two wintering seasons (from November to February) by faecal DNA metabarcoding using chloroplast trnL and mitochondrial CO1 region sequences. We examined 420 faecal samples and found different patterns of dietary diversity and composition among the duck species. The pattern also differed between plant and invertebrate food. Dietary niche partitioning was clear in plant food. Large-bodied ducks intensively use crop plants, and other ducks might mediate competition by using terrestrial and aquatic plants that are suitable for their foraging behaviours or microhabitats. Dietary segregation among species was the most apparent in February, when the abundance of foraging ducks was the largest. This study illustrated the complex pattern of dietary niche partitioning of dabbling ducks in agricultural landscapes, which might be difficult to evaluate by conventional approaches. The availability of crop plants, as well as other plant food resources in flooted areas and farmland dikes, may enable ducks to coexist by spatial or behavioural resource partitioning.

Genetic and demographic signatures accompanying the evolution of the selfing syndrome in Daphne kiusiana, an evergreen shrub.

BACKGROUND AND AIMS: The evolution of mating systems from outcrossing to self-fertilization is a common transition in flowering plants. This shift is often associated with the 'selfing syndrome', which is characterized by less visible flowers with functional changes to control outcrossing. In most cases, the evolutionary history and demographic dynamics underlying the evolution of the selfing syndrome remain poorly understood. METHODS: Here, we characterize differences in the demographic genetic consequences and associated floral-specific traits between two distinct geographical groups of a wild shrub, Daphne kiusiana, endemic to East Asia; plants in the eastern region (southeastern Korea and Kyushu, Japan) exhibit smaller and fewer flowers compared to those of plants in the western region (southwestern Korea). Genetic analyses were conducted using nuclear microsatellites and chloroplast DNA (multiplexed phylogenetic marker sequencing) datasets. KEY RESULTS: A high selfing rate with significantly increased homozygosity characterized the eastern lineage, associated with lower levels of visibility and herkogamy in the floral traits. The two lineages harboured independent phylogeographical histories. In contrast to the western lineage, the eastern lineage showed a gradual reduction in the effective population size with no signs of a severe bottleneck despite its extreme range contraction during the last glacial period. CONCLUSIONS: Our results suggest that the selfing-associated morphological changes in D. kiusiana are of relatively old origin (at least 100 000 years ago) and were driven by directional selection for efficient self-pollination. We provide evidence that the evolution of the selfing syndrome in D. kiusiana is not strongly associated with a severe population bottleneck.

Effects of fruit dimorphism on genetic structure and gene flow in the coastal shrub Scaevola taccada.

BACKGROUND AND AIMS: Plant propagules often possess specialized morphologies that facilitate dispersal across specific landscapes. In the fruit dimorphism of a coastal shrub, Scaevola taccada, individual plants produce either cork-morph or pulp-morph fruits. The former is buoyant and common on sandy beaches, whereas the latter does not float, is bird-dispersed, and is common on elevated sites such as slopes on sea cliffs and behind rocky shores. We hypothesized that beach populations bridge the heterogeneous landscapes by serving as a source of both fruit types, while dispersal is biased for the pulp morph on elevated sites within the islands and for the cork morph between beaches of different islands. Based on this hypothesis, we predicted that populations in elevated sites would diverge genetically over time due to isolation by distance, whereas beach populations would maintain high genetic similarity via current gene flow. METHODS: The genetic structure and gene flow in S. taccada were evaluated by investigating genome-wide single nucleotide polymorphisms in plants from 17 sampling sites on six islands (belonging to the Ryukyu, Daito and Ogasawara Islands) in Japan. KEY RESULTS: Geographical isolation was detected among the three distant island groups. Analyses within the Ryukyu Islands suggested that sandy beach populations were characterized by genetic admixture, whereas populations in elevated sites were relatively isolated between the islands. Pairwise FST values between islands were lowest between sandy beaches, intermediate between sandy beaches and elevated sites, and highest between elevated sites. CONCLUSIONS: Dispersal across the ocean by cork morphs is sufficiently frequent to prevent genetic divergence between beaches of different islands. Stronger genetic isolation of elevated sites between islands suggests that bird dispersal by pulp morphs is restricted mainly within islands. These contrasting patterns of gene flow realized by fruit dimorphism provide evidence that fruit characteristics can strongly mediate genetic structure.

Demography and selection analysis of the incipient adaptive radiation of a Hawaiian woody species.

Ecological divergence in a species provides a valuable opportunity to study the early stages of speciation. We focused on Metrosideros polymorpha, a unique example of the incipient radiation of woody species, to examine how an ecological divergence continues in the face of gene flow. We analyzed the whole genomes of 70 plants collected throughout the island of Hawaii, which is the youngest island with the highest altitude in the archipelago and encompasses a wide range of environments. The continuous M. polymorpha forest stands on the island of Hawaii were differentiated into three genetic clusters, each of which grows in a distinctive environment and includes substantial genetic and phenotypic diversity. The three genetic clusters showed signatures of selection in genomic regions encompassing genes relevant to environmental adaptations, including genes associated with light utilization, oxidative stress, and leaf senescence, which are likely associated with the ecological differentiation of the species. Our demographic modeling suggested that the glaberrima cluster in wet environments maintained a relatively large population size and two clusters split: polymorpha in the subalpine zone and incana in dry and hot conditions. This ecological divergence possibly began before the species colonized the island of Hawaii. Interestingly, the three clusters recovered genetic connectivity coincidentally with a recent population bottleneck, in line with the weak reproductive isolation observed in the species. This study highlights that the degree of genetic differentiation between ecologically-diverged populations can vary depending on the strength of natural selection in the very early phases of speciation.

The genome of Shorea leprosula (Dipterocarpaceae) highlights the ecological relevance of drought in aseasonal tropical rainforests.

Hyperdiverse tropical rainforests, such as the aseasonal forests in Southeast Asia, are supported by high annual rainfall. Its canopy is dominated by the species-rich tree family of Dipterocarpaceae (Asian dipterocarps), which has both ecological (e.g., supports flora and fauna) and economical (e.g., timber production) importance. Recent ecological studies suggested that rare irregular drought events may be an environmental stress and signal for the tropical trees. We assembled the genome of a widespread but near threatened dipterocarp, Shorea leprosula, and analyzed the transcriptome sequences of ten dipterocarp species representing seven genera. Comparative genomic and molecular dating analyses suggested a whole-genome duplication close to the Cretaceous-Paleogene extinction event followed by the diversification of major dipterocarp lineages (i.e. Dipterocarpoideae). Interestingly, the retained duplicated genes were enriched for genes upregulated by no-irrigation treatment. These findings provide molecular support for the relevance of drought for tropical trees despite the lack of an annual dry season.

Development of microsatellite markers for the geographically parthenogenetic stick insect Phraortes elongatus (Insecta: Phasmatodea).

Many plant and animal species exhibit geographic parthenogenesis, wherein unisexual (= parthenogenetic) lineages are more common in their marginal habitats such as high latitude or altitudes than their closely related bisexual counterparts. The Japanese stick insect, Phraortes elongatus (Thunberg) (Insecta: Phasmatodea), is known as a geographically parthenogenetic species due to the existence of both bisexual and unisexual populations. Here, we developed microsatellite markers to infer the genetic variation among populations of P. elongatus. Totally, 13 primer pairs were developed for the species, and they were tested on 47 samples collected from both a bisexual population and a unisexual population. All 13 loci were polymorphic in the bisexual population, whereas no loci were polymorphic in the unisexual population. The loss of variation in the unisexual population implies automixis with terminal fusion or gamete duplication as the mode of parthenogenesis. The markers developed in this study will be helpful for further comprehensive analysis of the genetic diversity and gene flow between bisexual and parthenogenetic lineages of P. elongatus.

Directional seed and pollen dispersal and their separate effects on anisotropy of fine-scale spatial genetic structure among seedlings in a dioecious, wind-pollinated, and wind-dispersed tree species, Cercidiphyllum japonicum.

Prevailing directions of seed and pollen dispersal may induce anisotropy of the fine-scale spatial genetic structure (FSGS), particularly in wind-dispersed and wind-pollinated species. To examine the separate effects of directional seed and pollen dispersal on FSGS, we conducted a population genetics study for a dioecious, wind-pollinated, and wind-dispersed tree species, Cercidiphyllum japonicum Sieb. et Zucc, based on genotypes at five microsatellite loci of 281 adults of a population distributed over a ca. 80 ha along a stream and 755 current-year seedlings. A neighborhood model approach with exponential-power-von Mises functions indicated shorter seed dispersal (mean = 69.1 m) and much longer pollen dispersal (mean = 870.6 m), effects of dispersal directions on the frequencies of seed and pollen dispersal, and the directions with most frequent seed and pollen dispersal (prevailing directions). Furthermore, the distance of effective seed dispersal within the population was estimated to depend on the dispersal direction and be longest at the direction near the prevailing direction. Therefore, patterns of seed and pollen dispersal may be affected by effective wind directions during the period of respective dispersals. Isotropic FSGS and spatial sibling structure analyses indicated a significant FSGS among the seedlings generated by the limited seed dispersal, but anisotropic analysis for the seedlings indicated that the strength of the FSGS varied with directions between individuals and was weakest at a direction near the directions of the most frequent and longest seed dispersal but far from the prevailing direction of pollen dispersal. These results suggest that frequent and long-distance seed dispersal around the prevailing direction weakens the FSGS around the prevailing direction. Therefore, spatially limited but directional seed dispersal would determine the existence and direction of FSGS among the seedlings.

Inferring historical survivals of climate relicts: the effects of climate changes, geography, and population-specific factors on herbaceous hydrangeas.

Climate relicts hold considerable importance because they have resulted from numerous historical changes. However, there are major interspecific variations among the ways by which they survived climate changes. Therefore, investigating the factors and timing that affected population demographics can expand our understanding of how climate relicts responded to historical environmental changes. Here, we examined herbaceous hydrangeas of genus Deinanthe in East Asia, which show limited distributions and a remarkable disjunction between Japan and central China. Chloroplast genome and restriction site-associated DNA sequencing revealed that speciation event occurred in the late Miocene (ca. 7-9 Mya) in response to global climate change. Two lineages apparently remained not branched until the middle Quaternary, and afterwards started to diverge to regional population groups. The narrow endemic species in central China showed lower genetic diversity (He = 0.082), as its population size rapidly decreased during the Holocene due to isolation in montane refugia. Insular populations in the three Japanese islands (He = 0.137-0.160) showed a genetic structure that was inconsistent with sea barriers, indicating that it was shaped in the glacial period when its range retreated to coastal refugia on the exposed sea floor. Demographic modelling by stairway-plot analysis reconstructed variable responses of Japanese populations: some experienced glacial bottlenecks in refugial isolation, while post-glacial range expansion seemingly exerted founder effects on other populations. Overall, this study demonstrated the involvement of not just one, but multiple factors, such as the interplay between climate changes, geography, and other population-specific factors, that determine the demographics of climate relicts.

A Narrow Endemic or a Species Showing Disjunct Distribution? Studies on Meehania montis-koyae Ohwi (Lamiaceae).

Meehania montis-koyae Ohwi (Lamiaceae), which has been considered a narrow endemic and endangered species in Japan, was found in eastern China in 2011. China and Japan belong to the same floristic region and share many plant species, but it is very rare that Japanese narrow endemic species are newly found outside of the country. We examined herbarium specimens of both countries, and conducted analyses of molecular phylogenetics, population genetics, and divergence time estimation using two nuclear (ITS and ETS) gene regions and MIG-seq data. Chinese plants tend to become larger than Japanese, and they are different in leaf shape and floral features. Molecular phylogenetic analysis shows Chinese and Japanese M. montis-koyae are the closest relatives to each other. Population genetic analysis indicates no current gene flow between the Chinese and Japanese populations, and divergence time analysis shows they were separated during the late Miocene. We reach the conclusion that Chinese and Japanese M. montis-koyae have already become distinct biological entities, and a new taxon name Meehania zheminensis A. Takano, Pan Li, G.-H.Xia is proposed for the Chinese plants. A key to Asian Meehania species is provided.

Genomic insights on the contribution of balancing selection and local adaptation to the long-term survival of a widespread living fossil tree, Cercidiphyllum japonicum.

'Living fossils' are testimonies of long-term sustained ecological success, but how demographic history and natural selection contributed to their survival, resilience, and persistence in the face of Quaternary climate fluctuations remains unclear. To better understand the interplay between demographic history and selection in shaping genomic diversity and evolution of such organisms, we assembled the whole genome of Cercidiphyllum japonicum, a widespread East Asian Tertiary relict tree, and resequenced 99 individuals of C. japonicum and its sister species, Cercidiphyllum magnificum (Central Japan). We dated this speciation event to the mid-Miocene, and the intraspecific lineage divergence of C. japonicum (China vs Japan) to the Early Pliocene. Throughout climatic upheavals of the late Tertiary/Quaternary, population bottlenecks greatly reduced the genetic diversity of C. japonicum. However, this polymorphism loss was likely counteracted by, first, long-term balancing selection at multiple chromosomal and heterozygous gene regions, potentially reflecting overdominance, and, second, selective sweeps at stress response and growth-related genes likely involved in local adaptation. Our findings contribute to a better understanding of how living fossils have survived climatic upheaval and maintained an extensive geographic range; that is, both types of selection could be major factors contributing to the species' survival, resilience, and persistence.

Corrigendum: Development of microsatellite markers for the annual andromonoecious herb Commelina communis f. ciliata (Commelinaceae) [Genes Genet. Syst. (2019) 94, p. 133-138].

Table 1 on p. 135 should be replaced with the corrected Table 1 shown bellow.

Corrigendum: New microsatellite markers recognize differences in tandem repeats among four related Gastrodia species (Orchidaceae) [Genes Genet. Syst. (2019) 94, p. 225-229].

Tables 3 and 4 on p. 228 should be replaced with the corrected Tables 3 and 4 shown bellow.

Genetic consequences of plant edaphic specialization to solfatara fields: Phylogenetic and population genetic analysis of Carex angustisquama (Cyperaceae).

Edaphic specialization is one of the main drivers of plant diversification and has multifaceted effects on population dynamics. Carex angustisquama is a sedge plant growing only on heavily acidified soil in solfatara fields, where only extremophytes can survive. Because of the lack of closely related species in similar habitats and its disjunct distribution, the species offers ideal settings to investigate the effects of adaptation to solfatara fields and of historical biogeography on the genetic consequences of plant edaphic specialization to solfatara fields. Here, genome-wide single nucleotide polymorphisms were used to reveal the phylogenetic origin of C. angustisquama, and 16 expressed sequence tag-simple sequence repeat markers were employed to infer population demography of C angustisquama. Molecular phylogenetic analysis strongly indicated that C. angustisquama formed a monophyletic clade with Carex doenitzii, a species growing on nonacidified soil in the sympatric subalpine zone. The result of population genetic analysis showed that C. angustisquama has much lower genetic diversity than the sister species, and notably, all 16 loci were completely homozygous in most individuals of C. angustisquama. Approximate Bayesian computation analysis supported the model that assumed hierarchical declines of population size through its evolutionary sequence. We propose that the edaphic specialist in solfatara fields has newly attained the adaptation to solfatara fields in the process of speciation. Furthermore, we found evidence of a drastic reduction in genetic diversity in C. angustisquama, suggesting that the repeated founder effects associated with edaphic specialization and subsequent population demography lead to the loss of genetic diversity of this extremophyte in solfatara fields.

New microsatellite markers recognize differences in tandem repeats among four related Gastrodia species (Orchidaceae).

Gastrodia is the most species-rich genus among mycoheterotrophic plants, and is thus an essential taxon to understand the mechanism of species diversification in mycoheterotrophs. In this study, we developed microsatellite markers with high transferability for four Gastrodia species to examine genetic differentiation and similarity among species, populations and individuals. The 12 microsatellite markers developed from a G. fontinalis library showed high transferability for the ramets that identified G. nipponica, G. kuroshimensis and G. takeshimensis. In addition to the high transferability of these markers, we observed low allele variation within a sampled population of each species and allele differences among the four species. The 12 markers described here will be useful for investigating the genetic differences among and within the Gastrodia species, which evolved by a limitation of gene flow.

Development of microsatellite markers for the endangered sleeper Eleotris oxycephala (Perciformes: Eleotridae).

The amphidromous sleeper Eleotris oxycephala (Perciformes: Eleotridae) is mainly distributed along the Kuroshio Current in East Asia, and this current is thought to be the main driver of the species' dispersal. Due to anthropogenic environmental changes in rivers, E. oxycephala is ranked as a threatened or near-threatened species in the red lists of 12 prefectures in Japan. Moreover, there is concern that the species' dispersal pattern could be changed due to fluctuations in the Kuroshio Current caused by global warming. In this study, 40 microsatellite markers were developed for E. oxycephala, and their suitability was tested on 43 individuals from two populations of E. oxycephala from Kanagawa and Miyazaki Prefectures. The number of alleles, expected heterozygosity and fixation index at each locus were 2-10 (mean = 5.350), 0.034-0.860 (mean = 0.650) and -0.261-0.448 (mean = 0.065), respectively. Furthermore, there was a lack of genetic difference between the two populations (FST = 0.008, F'ST = 0.024), indicating widespread gene flow via the Kuroshio Current. These markers will be useful to evaluate the genetic structure and infer population demographic history of E. oxycephala populations, which may assist in the conservation of this species.

Development of microsatellite markers for the annual andromonoecious herb Commelina communis f. ciliata (Commelinaceae).

Commelina communis f. ciliata (Commelinaceae), a newly distinguished taxon, is an annual andromonoecious herb exhibiting a mixed mating system, the details of which remain unclear. We developed microsatellite markers for use in exploring the evolution of andromonoecy and mixed mating in the species. Fifteen microsatellite loci were developed using next-generation sequencing. The primer sets were used to evaluate 65 C. communis f. ciliata individuals from three populations in Japan; we found 1-13 alleles per locus and the expected heterozygosity ranged from 0.00 to 0.76. The markers are potentially useful to examine intra- and interspecies genetic structure and the mixed mating strategy of Commelina species via paternity analysis.

Endangered island endemic plants have vulnerable genomes.

Loss of genetic diversity is known to decrease the fitness of species and is a critical factor that increases extinction risk. However, there is little evidence for higher vulnerability and extinction risk in endangered species based on genomic differences between endangered and non-endangered species. This is true even in the case of functional loci, which are more likely to relate to the fitness of species than neutral loci. Here, we compared the genome-wide genetic diversity, proportion of duplicated genes (PD), and accumulation of deleterious variations of endangered island endemic (EIE) plants from four genera with those of their non-endangered (NE) widespread congeners. We focused on exhaustive sequences of expressed genes obtained by RNA sequencing. Most EIE species exhibited significantly lower genetic diversity and PD than NE species. Additionally, all endangered species accumulated deleterious variations. Our findings provide new insights into the genomic traits of EIE species.

Environmental DNA Collected from Snow Tracks is Useful for Identification of Mammalian Species.

Noninvasive genetic analysis is being used increasingly in field surveys. However, detecting large and middle-sized mammals, such as Carnivora species, using noninvasive samples, such as scat or hair, is time- and labor-intensive due to their low densities and elusive behaviors. As snow tracks are the most frequently encountered natural signs of terrestrial mammals in winter, we employed several methods to recover environmental DNA (eDNA) from snow tracks. We performed both DNA metabarcoding and Sanger sequence analyses, in combination with universal primers on the mitochondrial 12S rRNA gene for mammals and taxon-specific primers on the mitochondrial NADH dehydrogenase subunit 2 gene for Martes species (martens and sables in Mustelidae). Snow samples of four Martes melampus tracks, one Cervus nippon track, one Vulpes vulpes track, and the track of an unidentified Carnivora species were collected from a snowfall area in Kyoto, Japan, in February 2018. Regarding DNA metabarcoding analyses, the sequences of three Carnivora species (M. melampus, V. vulpes, and Canis lupus familiaris) and a deer (C. nippon) were obtained from their respective snow tracks. Using Sanger sequencing, eDNA on snow tracks was recovered at the species level except for M. melampus using universal primers, while eDNA of M. melampus was sequenced using Martes-specific primers. Snow track surveys in combination with eDNA techniques could dramatically improve the efficiency of monitoring and conservation of mammals.